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Determination of the Code Phase Between a 
Code Modulated Signal and a Replica Code Sequence 

CROSS REFERENCE TO RELATED APPLICATIONS 

5 

This application claims priority under 35 USC §119 to 
International Patent Application No. PCT/IB02/04421 filed 
on October 24, 2002. 

10 FIELD OF THE INVENTION 

The invention relates to a method for determining the code 
phase between a code modulated signal received at a 
receiver and an available replica code sequence. The 
15 invention relates equally to a receiver, to an electronic 

device and a communication system comprising a receiver and 
to a device communicating with a receiver. 

BACKGROUND OF THE INVENTION 

20 

The code phase between a code modulated signal received at 
a receiver and an available replica code sequence has to be 
determined for example for CDMA (Code Division Multiple 
Access) spread spectrum receivers. 

25 

For a spread spectrum communication in its basic form, a 
data sequence is used by a transmitting unit to modulate a 
sinusoidal carrier and then the bandwidth of the resulting 
signal is spread to a much larger value. For spreading the 

30 bandwidth, the single -frequency carrier can be multiplied 

for example by a high-rate binary pseudo-random noise (PRN) 
code sequence comprising values of -1 and 1, which code 
sequence is known to a receiver. Thus, the signal that is 
transmitted includes a data component, a PRN component, and 

35 a sinusoidal carrier component. A PRN code period comprises 
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typically 1023 chips, the term chips being used to 
designate the bits of the code conveyed by the transmitted 
signal, as opposed to the bits of the data sequence. 

5 A well known system which is based on the evaluation of 
such code modulated signals is GPS (Global Positioning 
System) . In GPS, code modulated signals are transmitted by 
several satellites that orbit the earth and received by GPS 
receivers of which the current position is to be 

10 determined. Each of the satellites transmits two microwave 
carrier signals. One of these carrier signals LI is 
employed for carrying a navigation message and code signals 
of a standard positioning service (SPS) . The LI carrier 
signal is modulated by each satellite with a different C/A 

15 (Coarse Acquisition) Code known at the receivers. Thus, 

different channels are obtained for the transmission by the 
different satellites. The C/A code, which is spreading the 
spectrum over a 1 MHz bandwidth, is repeated every 1023 
chips, the epoch of the code being 1 ms . The carrier 

20 frequency of the LI signal is further modulated with the 
navigation information at a bit rate of 50 bit/s. The 
navigation information, which constitutes a data sequence, 
can be evaluated for example for determining the position 
of the respective receiver. 

25 

A receiver has to have access to a synchronized replica of 
the modulation code which was employed for a received code 
modulated signal, in order to be able to de- spread the data 
sequence of the signal. To this end, a synchronization has 

30 to be performed between the received code modulated signal 
and an available replica code sequence. Usually, an initial 
synchronization called acquisition is followed by a fine 
synchronization called tracking. In both synchronization 
scenarios, a correlator is used to find the best match 

35 between the replica code sequence and the received signal 
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and thus to find their relative shift called code phase. 
The search can be performed with different assumptions on 
an additional frequency modulation of the received signal. 
Such an additional modulation may occur for example due to 
5 a Doppler effect and/or a receiver clock inaccuracy and can 
be as large as +/-6 kHz. 

Two main types of correlators have been suggested so far. A 
first type of correlators performs a direct correlation of 

10 the replica code sequence and the received signal in the 

time domain. This implies that a dedicated processing step 
is carried out for each possible code phase. In case there 
is a large number of code phases to check, the 
computational burden is significant, especially for 

15 software based receivers. There exist different 
implementation approaches for the first type of 
correlators, which may be formed with matched filters or 
ordinary correlators. A second type of correlator relies on 
frequency domain acquisition techniques employing e.g. 

20 Discrete Fourier Transforms (DFT) , which enable a parallel 
processing for all possible code phases and thus a 
reduction of the computational burden. 

Figure 1 illustrates a known DFT based circular correlation 
25 in the frequency domain. To simplify the illustration, the 
modulation code is supposed to comprise eight samples. In 
practice, the code will usually comprise a larger number of 
samples, e.g. 1024 samples. First, a vector 11 with eight 
samples of a received code modulated signal is provided to 
30 the correlator. Each sample in figure 1 is indicated by a 
small circle. The correlator performs a DFT 12 of the 
provided vector 11, resulting in another vector 13 with 
eight samples. Further, the correlator retrieves or 
calculates a conjugate 14 of the DFT of a vector comprising 
35 eight samples of an available replica code sequence. The 



Express Mail No. EV30371 1515US 3 



PATENT 

Attorney Docket No. 915-007.051 

DFT vector 13 of the received signal and the conjugate 14 
of the DFT vector of the replica code sequence are then 
multiplied pointwise 15. For the resulting vector 16 of 
again eight samples, an Inverse Discrete Fourier Transform 
5 (IDFT) 17 is performed, which results again in a vector 18 
comprising eight samples. Each sample of the output IDFT 
vector 18 corresponds to a correlation value for another 
one of all possible circular shifts. The vector may 
comprise for example the sample values [0.5 7.8 2.3 5.3 2.9 

10 3.4 4.5 0.7] which are associated in this order to the code 
phases [01234567], In the presented example, the 
maximal value of the output samples is 7.8, thus the found 
code phase is 1. This means that the replica code is 
shifted by one sample relative to the received code of the 

15 code modulated signal. 

In principle, the phase of the received code relative to 
the available replica code sequence can have any possible 
value. In some situations, however, the range of the 

20 possible code phases can be reduced based on some apriori 
knowledge regarding e.g. the position of the transmitting 
unit, the position of the receiver and the time of 
transmission of the received signal. Such apriori knowledge 
may be available for example at assisted GPS receivers (A- 

25 GPS) . Assisted GPS receivers use additional information, 
provided e.g. by a cellular network, to accelerate and 
simplify the algorithms used for position calculations. 

When the location of a receiver is already known with a 
30 certain accuracy in addition to available ephemeris and 
time information, the synchronization procedures for 
acquisition and tracking would advantageously not check all 
possible values of the code phases but only a limited 
number. The conventional search in GPS is carried out for 
35 1024 chips, which corresponds to an uncertainty area of 
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around 3 00km. Certain scenarios on the newly designed 
Galileo system, the European analog of GPS, could even have 
a search uncertainty area of a few thousands of kilometers. 
In an urban area, though, the position of the receiver 
5 might be known with an accuracy of about 1km, e.g. from 
some assistance. This knowledge may be exploited for 
performing only a limited search. 

There are several situations in which the range of possible 

10 code phases can be limited. For example, if a specific GPS 
satellite is acquired and tracked and the position of the 
GPS receiver is known with an accuracy of about 5 0km, then 
the phase uncertainty is limited to l/6 th of the whole range 
of 1023 possible code phases, as the GPS time can be 

15 reconstructed with a good accuracy. Further, if a GPS 
satellite was tracked and the position of the receiver 
determined, and then the signal is lost again, the GPS time 
will still continue to be quite accurate, since the 
internal clock was recently initialized accurately. In an 

20 urban area, it can further be assumed that the speed of the 
receiver is limited to 50km/h, i.e. to about 20m/s. Thus, 
the receiver can be assumed to be in a 20km area from the 
previously determined position for around 20min, and the 
phase uncertainty is limited l/10 th of the whole range of 

25 1023 code phases. In the latter case, the receiver might 

even know without assistance that only a limited number of 
code phases is possible. 

Currently, however, a limited search of code phases can 
30 only be realized with correlators performing a correlation 
in the time domain. Known DFT based methods inherently 
perform the search of all possible code phase in parallel. 
Therefore, their usage is not feasible in situations in 
which the search is to be carried out only over a limited 
35 number of all possible code phases. With conventional DFT 
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correlators, known limitations for the code phase can only 
be evaluated after the IDFT. Thus, it is a disadvantage of 
conventional DFT correlators that they perform in many 
situations unnecessary computations. Depending on the 
5 extent to which the range of the possible code phases can 
be limited, the use of correlators operating in the time 
domain might even be more reasonable again. 

SUMMARY OF THE INVENTION 

It is an object of the invention to reduce the amount of 
required processing in a time to frequency transform based 
correlation procedure, which is employed for determining 
the code phase between a received code modulated signal and 
an available replica code sequence. It is in particular an 
object of the invention to reduce the amount of required 
processing for the case that the number of possible code 
phases can be restricted beforehand. 

20 A method is proposed which comprises as a first step 

performing a multiplication between samples of a first 
vector and samples of a second vector resulting in a third 
vector. This multiplication can be realized for instance as 
elementwise or pointwise multiplication. The first vector 

25 is generated based on the received code modulated signal in 
an operation including a time to frequency transform, and 
the second vector is generated based on the replica code 
sequence in an operation including a time to frequency 
transform. It is to be noted that the actual generation of 

30 the second vector does not necessarily constitute a part of 
the proposed method. It can be stored for example for each 
available replica code sequence. Then, the obtained third 
vector is divided into sections, and the samples in each 
section are summed. Out of the summed samples, a reduced 

35 fourth vector is formed. Finally, a frequency to time 
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transform of said fourth vector is performed. The frequency 
to time transform results in a fifth vector. Each sample of 
this fifth vector represents a correlation value for a 
different code phase between the received code modulated 
5 signal and the available replica code sequence. 

Moreover, a receiver, an electronic device comprising a 
receiver and some other device are proposed, either 
comprising means for carrying out the steps of the proposed 
method. In case the processing is performed in another unit 
than the receiver, the required information about the 
received signals is forwarded by the receiver to this unit. 
The proposed other device can be for instance a network 
element of a network. The object is also reached with a 
system comprising a receiver and a device, in which system 
either the receiver or the device comprises means for 
carrying out the steps of the proposed method. In case the 
receiver performs the processing, the device may provide 
assistance data to the receiver. 

The invention proceeds from the idea that the calculations 
performed for those code phases that do not lie within a 
limited range of possible code phases do not have to be 
skipped only in the frequency to time transform itself. 
Instead, the vector for which the frequency to time 
transform is determined can advantageously be reduced 
beforehand . 

A time to frequency transform has the useful property that 
30 a circular shift in the input vector results in a complex 
sinusoidal modulation of the transform outputs which are 
obtained in case there is no shift. Thus, the transform 
outputs are the same for all possible shifts, except that 
they are modulated differently. The modulation frequency 
35 depends on the shifting distance, i.e. the larger the 
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shifting, the higher the modulation frequency. In case, for 
example, the outputs of a time to frequency transform of a 
received signal are multiplied with the output of a time to 
frequency transform of an inverted conjugate of the replica 
5 code sequence, components of the correlation in the 
frequency domain modulated according to the shift are 
obtained. The subsequent frequency to time transform 
detects this modulation and outputs the largest value at a 
vector index corresponding to the shift value. If now the 

10 range of the possible code phases is restricted to a known 
value, this means that the modulation in the frequency 
domain is also restricted. Thus, it is possible to 
integrate correlation components already in the frequency 
domain without a preceding demodulation by a frequency to 

15 time transform. As a result, only values for those code 

phases which are closest to the alignment are output . The 
correct code phase is the output index which has the 
largest output value. The integration length should depend 
on the range of possible code phases and defines the 

20 modulation frequency range in the frequency domain. 

The invention thus modifies the known time to frequency 
transform based correlation method to allow a parallel 
search over a restricted range of possible code phases. 
25 With the proposed modification, the size and complexity of 
the frequency to time transform can be reduced in certain 
scenarios, which enables an optimization of the frequency 
domain computations . 

30 Compared with conventional time to frequency transform 
based methods, the complexity may be reduced in some 
situations up to tens or even hundreds of times. In case of 
1023 possible code phases, for instance, the conventional 
time to frequency transform based frequency domain 

35 technique searches over all 1023 possibilities, while the 
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invention is suited to optimize the frequency domain 
processing by reducing the search to e.g. 16 or 32 code 
phases . 

5 The main complexity of a time to frequency transform based 
correlator is distributed equally between the forward and 
inverse transforms, and if the frequency to time transform 
size reduces dramatically, then the entire complexity will 
be reduced down to half. Moreover, in certain time to 

10 frequency transform based correlation methods it is 
possible to calculate the forward time to frequency 
transform only one time and to use the result with 
different replicas for different satellites and for 
different frequency bands by circularly shifting the 

15 replica. Such a method was proposed in D. Akopian, I. 

Kontola, H. Valio, S. Turunen, "Method in a receiver and a 
receiver, 11 patent application, Nokia Mobile Phones, 1999, 
and by D. Akopian in "A fast satellite acquisition method", 
ION-GPS f 2001 Conference, Salt Lake City, USA, Sept. 11-14, 

20 2001. In this case, the main complexity of the time to 

frequency transform based correlator over many frequency 
and satellite searches is concentrated on the frequency to 
time transform stage. Therefore, reducing the frequency to 
time transform size dramatically, i.e. by orders of 

25 magnitude, will reduce the overall cost of the correlation 
stage by the same amount . 

The computational complexity reduction can be utilized by 
using a slower processor, resulting in a reduced power 
30 consumption or enabling a software-only implementation. On 
the other hand, with the same computational power it will 
be possible to perform algorithms with low complexity 
faster and thus to reduce delays . 
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Preferred embodiments of the invention become apparent from 
the dependent claims. 

The number of sections, into which the vector resulting in 
5 the multiplication is divided, is preferably selected based 
on an available information on a limited range of possible 
code phases. The number of sections should be equal to or 
larger than the number of possible code phases in this 
limited range. The limited range of possible code phases 
10 can be determined in particular based on available 
information on a position of the receiver. 

Advantageously, but not necessarily, the sections are of 
equal size. If they are not of equal size, the outputs will 
15 be distorted, but the frequency to time transform can be 
modified to account for this inequality. 

In order to cope with a multipath propagation of the code 
modulated signal, the code modulated signal may be 
20 correlated in accordance with the invention with a 

plurality of identical replica code sequences which are 
shifted in phase. To this end, a plurality of similar 
correlators may be provided. 

25 The first and the second vector multiplied in the 
multiplication can be obtained in various ways. 

The first vector can be obtained for example by performing 
a time to frequency transform of the received code 

30 modulated signal. In this case, the second vector can be 
given e.g. by a vector resulting in a time to frequency 
transform of the inverted conjugate of the replica code 
sequence. Alternatively, the second vector can be given in 
this case by the conjugate of a vector resulting in a time 

35 to frequency transform of the replica code sequence. 
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On the other hand, the second vector can be obtained by 
performing a time to frequency transform of the replica 
code sequence. In this case, the first vector can be given 
5 e.g. by a vector resulting in a time to frequency transform 
of the inverted conjugate of the received code modulated 
signal. Alternatively, the first vector can be given in 
this case by the conjugate of a vector resulting in a time 
to frequency transform of the received code modulated 
10 signal. 

The time to frequency transform performed for obtaining the 
first and second vector can be in particular, though not 
exclusively, a DFT . Correspondingly, the frequency to time 
15 transform performed for obtaining the fifth vector can be 
in particular, though not exclusively, an IDFT . 

The time to frequency transform can be implemented as a 
fast computational method, for example a Fast Fourier 
20 Transform or any other suitable approach. 

The invention can be used in both acquisition and tracking 
schemes. In tracking, e.g. multiple shifted correlators 
could be utilized for multipath mitigation. 

25 

The invention may be used in both cases for determining the 
code phase and the frequency of a remaining complex 
sinusoidal modulation, i.e. of the sinusoidal modulation 
which remains after the carrier has been wiped off from the 

30 received signal based on the known nominal carrier 

frequency. The code phase is determined according to the 
peaks of a cross-correlation function, and the correlation 
is calculated at initial code wipe-off stages. The 
processing for weak signals requires additional coherent 

35 and non-coherent integrations. The invention can therefore 
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also be used as a building block for other methods 
implementing different scenarios of coherent and/or non- 
coherent processing for possible multiple frequency- 
candidates . 

The invention can be implemented in hardware or in 
software. In case the invention is employed as part of 
acquisition and tracking algorithms, the implementation 
corresponds advantageously to the implementation of these 
algorithms . 

The invention can be employed in particular, though not 
exclusively, for CDMA spread spectrum receivers, for 
instance for a receiver of a positioning system like GPS or 
15 Galileo. 

BRIEF DESCRIPTION OF THE FIGURES 

Other objects and features of the present invention will 
20 become apparent from the following detailed description 

considered in conjunction with the accompanying drawings, 
wherein 

Fig. 1 illustrates a DFT based correlation according to 
25 the state of the art; and 

Fig. 2 illustrates a DFT based correlation according to 
an embodiment of the invention. 

DETAILED DESCRIPTION OF THE INVENTION 

30 

Figure 1 has already been described above. 

Figure 2 illustrates an exemplary embodiment of the method 
according to the invention implemented in an A-GPS 
35 receiver. The receiver comprises a receiving unit for 
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receiving signals from different GPS satellites which are 
modulated with different C/A-codes, each comprising 1023 
chips per code period. Moreover, the GPS receiver comprises 
a tracking unit with a correlator employing DFT based 
5 frequency domain acquisition techniques for acquiring and 
tracking received satellite signals. The tracking unit has 
access to a replica code sequence for each of the GPS 
satellites. The GPS receiver is furthermore included in a 
mobile terminal of a communication system. A 
10 microcontroller unit (MCU) of the GPS receiver is able to 
store and evaluate assistance information received by this 
mobile terminal from a communication network or information 
available at the receiver. 

15 It will be assumed in the following that the number of 

samples N per code period of received satellite signals is 
equal to a power of two, i.e. in the current example 1024 
chips instead of 1023 chips. Other cases can be explored in 
a similar manner. 

20 

In a first step, a C/A code modulated input signal 
x = {x 0 , . . . , received at the receiver is provided to 

the correlator. In figure 2, the input signal x is 
represented for reasons of simplicity by a vector 21 
25 comprising eight samples instead of 1024 samples. As in 
figure 1, each sample is indicated by a small circle in 
figure 2 . 

In the correlator, a DFT 22 of the input signal 
30 x = {x 0 , . . . is performed. In case the DFT matrix is 

denoted as F, the resulting vector y 1 is given by 

y 1 = f n* • 
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Vector y 1 is represented in figure 2 by another vector 23 
comprising eight samples. 

Further, a vector r resulting in a DFT of the inverted 
5 conjugate of the replica code sequence is provided. In 
figure 2, this vector r is represented by yet another 
vector 24 comprising eight samples. 

In a next step 25, vector y 1 is multiplied pointwise with 
10 vector r. The resulting vector y 2 is given by 

y 2 — y 1 . * r 

The pointwise operation is denoted as Also vector y 2 

15 is represented in figure 2 by a vector 26 comprising eight 
samples . 

So far, the processing by the correlator corresponds to the 
known DFT based processing described with reference to 
20 figure 1. 

In contrast to the known processing, however, the vector y 2 
resulting in the pointwise multiplication 24 is not 
subjected immediately to an IDFT 27. Rather, it is first 
25 divided into K sections 29 of equal size. The value of K is 
set by the MCU to the number of possible code phases, which 
number is determined by the MCU based on available 
assistance data. 

30 In the example of figure 2, eight different code phases [0 
12 3 4 5 6 7] may exists in the whole. Due to available 
assistance data on a reference position of the receiver and 
of the satellite transmitting the received signal, it is 
known that currently at the most a code shift by 1 in one 

35 direction and by 2 in the opposite direction may occur. 
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This corresponds e.g. to the four possible code phases [0 1 
6 7] , since the code phases are circular. "0 XX corresponds 
to an exactly aligned input signal and replica code 
sequence, "1" corresponds to an alignment shifted by one 
5 sample in one direction, W 7 W corresponds to an alignment 
shifted by one sample in the opposite direction, and "6* 
corresponds to an alignment shifted by two samples in the 
opposite direction. The limitation to four code phases 
results in a value of K = 4. Thus, in figure 2 the eight 
10 samples are divided into four sections 29 of two samples 

each. In GPS, this would corresponds to 512 sections of two 
samples each. 

The samples in each section are then summed 30. The result 
15 is a vector y 3 of a reduced size K. Vector y 3 is 

represented in figure 2 by a vector 31 comprising four 
samples . 

The IDFT 27 is now applied to this reduced vector 31 
20 according to the following equation: 

z = f;V 

Vector z is represented in figure 2 by a vector 28 
25 comprising again four samples. 

The interpretation of the output index of the IDFT 2 7 and 
thus of the expected correlation peak index is as follows. 
The IDFT algorithm finds the code phases around the aligned 
30 position corresponding to a code phase of n 0 w . The possible 
values are {- K/2, . . . , - 1,0,1, - - • , K/2 - l} . The values of the 
first K/2 samples of vector z correspond to positive shifts 
{0,1, . . . ,K/2 — l} while the values of the next K/2 samples of 
vector z correspond to negative shifts {- K/2, . . . , - l} . In 
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the above example, the outputs of the IDFT are thus 
associated to the phases [0167]. Proceeding from the 
exemplary correlation values presented with reference to 
figure 1, the output vector 28 in figure 2 would be [0.5 
5 7.8 4.5 0.7]. The final result is the same as in figure 1, 
i.e. the sample with the maximum value is the second one, 
and thus the code phase is 1. This time, however, the 
output samples which are not needed due to an apriori 
knowledge of limitations for the possible code phases are 
10 not calculated at all. 

The vector z resulting in the IDFT 27 can further be used 
for an additional coherent and/or non-coherent processing 
which is performed for handling low strength signals in 
15 noise. 

Also in cases in which the correlation should be performed 
for different values of a remaining complex sinusoidal 
modulation after the carrier wipe-off, there is no need to 

20 perform a DFT transform for each possible modulation 

frequency. Instead, a shift of the transformed replica code 
sequence could be used. In this case, a reduction in 
complexity of the IDFT 2 7 according to the invention has a 
particular benefit. For example, when the A-GPS receiver 

25 position is known with an accuracy of 3 000m, this 

corresponds to a search area of approximately 10 code 
phases. For simplicity, the closest power of two, 16, is 
taken for the number of candidate code phases to search. 
The complexity reduction will then be approximately 

30 (1024*10) / (16*4) = 160 times over the conventional DFT 
based correlator. 

It is to be noted that the described embodiment constitutes 
only one of a variety of possible embodiments of the 
35 invention. 
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